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This paper addresses the problem of controlling a nonlinear plant with a slow actuator 
using singular perturbation method. For the known plant-actuator cascaded system the 
proposed scheme achieves tracking of a given reference model with considerably less control 
demand than would otherwise result when using conventional design techniques. This is 
the consequence of excluding the small parameter from the actuator dynamics via time 
scale separation. The resulting tracking error is within the order of this small parameter. 
For the unknown system the adaptive counterpart is developed based on the prediction 
model, which is driven towards the reference model by the control design. It is proven that 
the prediction model tracks the reference model with an error proportional to the small 
parameter, while the prediction error converges to zero. The resulting closed-loop system 
with all prediction models and adaptive laws remains stable. The benefits of the approach 
are demonstrated in simulation studies and compared to conventional control approaches. 


I. Introduction 


In many applications, the plant to be controlled has much slower dynamics than the actuator through 
which it is being controlled. That is, the cascaded plant-actuator system can be described as 

x(t) = f(x(t))+g(x(t))y{t) 

ez(t) = a(z(t)) + b(z(t))u(t) 

y{t) = c{z(t)), (1) 

where x(t) € M" is the plant’s state, z(t) € R m is the actuator’s state, y(t) £ 1 is the actuator’s output, 
u(t) € M is the actuator’s input, and e is a small constant. For these types of systems, the singular pertur- 
bations method has been successfully used for control design. 6 This approach requires some interconnection 
conditions to be satisfied, which can be imposed directly on the systems under the consideration (See for 
example 16 ) or in terms of Lyapunov functions and their derivatives (See for example 13 ). The developed 
controllers typically achieve local results (stability or asymptotic stability) of closed-loop systems (See for 
example Refs. 2, 6,8 just to mention a few of them). However, for some systems, global exponential stability 
can be achieved. 1 As it is shown in Ref., 5 when the exponential stability is achieved, the actuator’s dynamics 
can be neglected and the actuator’s output can be viewed as a control input for the plant. Therefore, the 
control design can follow any known scheme that is suitable for the given plant dynamics. 

In some applications, the actuator dynamics can be substantially slower than the plant’s dynamics. That 
is, the cascaded plant-actuator system is now described by the equations 


x{t) 

= f{x(t))+g{x(t))y{t) 


z(t) 

= a(z(t)) + b(z(t))u(t) 


y(t) 

= c(z(t)). 

(2) 
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This is the case when the specific actuator, which is designed to control the plant in ordinary scenarios, 
fails, and the substitute is used to rescue the emergency situation. For example, the aircraft engine can 
be used for the directional stability and control of the aircraft when the control surface deflection system 
fails. However, the engine dynamics is known to be much slower than the aircraft’s rotational dynamics with 
the rudder input. For these types of problems, the actuator dynamics cannot be ignored, and the actual 
control input has to be designed to generate the necessary actuator output. In general, the actual control 
input depends is inversely proportional to the small parameter, thus resulting in high gain control signal. 
Therefore, to generate the necessary torque or force, the actuator may require a control signal exceeding the 
physical limits. 

The goal of this paper is to show that the application of the singular perturbation method based time 
scale separation can avoid demanding high magnitude control signals, and to develop a control algorithm 
that guarantees the stable tracking of a given reference command in the presence of modeling uncertainties 
and external disturbances. 

The rest of the paper is organized as follows. First, we give preliminaries for the singular perturbations 
method, then we apply the method to the known plant’s and actuator’s dynamics. Afterwards, we present 
time scale separation based control design for the unknown system and provide stability analysis. Throughout 
the paper bold symbols are used for vectors, capital letters for matrices and small letters for scalars. 

II. Preliminaries 

Consider the singular perturbation model (1) with initial conditions £c(0) = Xq, z(o) = Zq. We assume 
that the functions f(x) and g(x) are continuously differentiable on some open connected set D x C M”, and 
the functions a(z) and b(z) are continuously differentiable on some open connected set D~ C R m . Let the 
control input have the form u(t) = ip(t,x,z), where the function ip(t,x,z) and its first partial derivatives 
with respect to (t,x,z) are continuous and bounded on any compact subset of xD x x D z . 

The model is said to be in standard form if the algebraic equation 

a(z) + b(z)ip(t, x, z) = 0 (3) 

has isolated roots for all (t, x) € [0, oo) x D x . 

Let z = h(t,x ) be an isolated root of (3). The system 

x{t) = f(x(t)) + g(x(t))h(t , x(t)) , a:(0) = x Q (4) 

is called a reduced system. Denote the solution of (4) by x{t). 

Introduce a change of variables £ = z — h(t, x) and t = et. The system 

= a (£ + h{t,x)) + b(£ + h(t,x))ip(t,x,£ + h(t,x )) , £(0) = z 0 - h(0, x 0 ) (5) 

where (t, x) € [0, oo) x D x are treated as fixed parameters, is called a boundary layer system. 

In this paper, for the stability analysis of the closed-loop systems, we will use the following version of 
Theorem 11.2 from Ref. 5 

Theorem 1 Consider the singular perturbation problem for system (1) with initial conditions a:(0) = xq, 
z(o) = Zq and control input u(t) = ip(t,x,z). Let z = h(t,x) be an isolated root of the algebraic equation 
(3). Assume that the following conditions are satisfied 

• The functions f(x) and g{x) are continuously differentiable on D x ; 

• The function a(z), b(z) and c(z) are continuously differentiable on D z ; 

• The function ip(t, x, z) ant its first partial derivatives with respect to its arguments are continuous and 
bounded on any compact subset of D x x D z ; 

• The origin is an exponentially stable equilibrium point of the reduced system (4), 7t x is the region of 
attraction; 

• The origin is an exponentially stable equilibrium point of the boundary-layer system (5), uniformly in 
(t,x), 1Z Z is the region of attraction. 



Then, for each compact set tt x C 1Z X and fl z C 1Z Z there is a positive constant £* such that for all t > 0, 
Xq G Cl x , Zq — h(0,xo) G and 0 < £ < £* , the singidar perturbation problem (1) has a unique solution 
x(t, e), z(t, e) on [0, oo), and 


x(t, e) — x(t) = 0(e) (6) 

holds uniformly for all t. > 0. Moreover, there exists a time instance T(e) such that 

z(t,e)-h(t,x(t))=0(e) (7) 


holds uniformly for all t, > T(e). 


III. Problem Formulation 

Consider a single-input single-output system 

x p (t) = A p x p (t) + b p lf P (x p (t)) + g P {x p (t))y a (t) + d P (t)\ 

y P (t) = cjx p (t), (8) 

with the initial conditions *(0) = x p o, where x p G R" is the plant’s state, f p {x p ) and g p (x p ) are continuously 
differentiable unknown functions representing the modeling uncertainties, and y a {t) G R is the control effort 
generated by a slow actuator. It is reasonable to assume that the actuator’s dynamics are stable, minimum 
phase and of known relative degree r, and are described by the system 

£~ 1 X a (t) = A a X a (t) +b a [wu(t) + fa(y a (t),P a ) +d a (t)} 

Va{t) = cjx a (t) (9) 

with initial conditions a: a (0) = x a Q. Here, x a G R m is the actuator’s state, A a is a Hurwitz matrix, the 
transfer function G 0 (s) = cf (si — A a )~ 1 b a has stable zeros, e is a small parameter, u> G (0, 1] represents 
the possible reduction in actuator’s control effectiveness, f a (ya,P a ) is a continuously differentiable function 
representing the modeling uncertainties, u(t) is the actuator’s control input, and time signals d p {t) and d a (t) 
represent unknown external disturbances. 

We notice that the actuator’s dynamics do not depend on the plant’s state to be controlled. However, 
the dynamics may depend on other external variables, which are lumped into the parameter p a . In general, 
the actuator’s state is not available for feedback. Therefore, the actuator related part of the overall control 
design must be considered in output feedback framework, assuming that the signal y a {t) is available for 
feedback. 

To guarantee the controllability of the plant, we also assume that the pair ( A p , b p ) is controllable, the 
pair ( Ap,c p ) is observable, and the function g p {x p ) never crosses zero. Without loss of generality, we take 
g p {x p ) as a positive function, that is, there exists a positive constant go such that g p (x p ) > go > o. 

The control objective is to design an input signal u(t) for the actuator such that the system’s output 
y p (t) tracks the output of a given reference model 

Xm(t) A m x m (t) + b m y c (t) 

y m (t) = c T x m (t), (10) 

with initial conditions x m (0) = x m o, where y c {t) is an external command that is bounded with its first 
time derivative. Here, A rn is a Hurwitz matrix of the form A, m = A p — b p k T for some control gain k , and 
b rn = vbp with v = (cj A~fb p ) 1 . Since A rn is Hurwitz, there exist a symmetric positive definite matrix P 
such that 


+ P -Am — ~Q 


( 11 ) 


for some symmetric positive definite matrix Q. 



IV. Control design for known systems 


In this section, we show that the magnitude of a control signal designed by the proposed time scale 
separation based algorithm is smaller than the magnitude of a control signal generated by the conventional 
algorithms by a term that is of the order of e _1 . For comparison, we use backstepping design, 7 which is well 
suited for the cascaded system comprised of the dynamics in (8) and the actuator dynamics (9). To this 
end we assume that the systems under consideration are known and disturbance free. Also, for the purposes 
of this section we assume that the functions f p (x p ), g p (x p ) and g a (ya) are sufficiently smooth, and the 
reference command and its sufficiently many time derivatives are bounded. In this case, for convenience we 
write the actuator dynamics in the following normal form 10 

£ _1 *a(0 = A) xl(t)+b 0 [-kjx a (t)+v(t)\ 

= A c 0 C(t) + Kx c a (t) 

y a (t ) = c 0 x c a (t), ( 12 ) 

where x c a £ R r is the portion of the actuator’s state corresponding to the chain of r integrators represented 
by the triplet (ej” £ R r , A G K rxr , bo € R r ) having the form 
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g jjm-r i s the actuator’s internal state, and the eigenvalues of Ag £ jg)m-r)x(m-r) are th e stable zeros 
of the actuator’s transfer function G a {s). For the notational simplicity, we introduce a new control variable 
v(t) = bum(t) + bg a {ya(t),Pa)i where b ^ 0 is the high frequency gain. Since the function g a (x a ,p a ) and 
parameters w and b are assumed to be known, u(t) can be readily computed after designing v(t). 

Let the tracking error be e(t) = x p (t) — x m (t). Its dynamics can be written as 


e{t) 


A m e(t) + b p k T x p (t ) + f p (x p (t )) + g p {x p (t))y a (t) - vy c (t) 


(13) 


Since for the error dynamics (13) there exists a continuously differentiable feedback control law in the form 


Va(t) = on(t, x p ) = g p 1 (x p (t)) [-k T x p (t) - fp(x p (t)) - vy c {t) 
that guarantees that the derivative of the radially unbounded Lyapunov function 

Vi (e(t)) = e T (t)Pe(t) 


(14) 


(15) 


is negative definite for all e £ R” 

Vi(£) = -e 1 (t)Qe(t) < 0, (16) 

the linear block backstepping lemma (Lemma 2.23) from Ref. 7 can be applied. To be prepared for that 
application, we introduce change of variables z c (t ) = T c x c a (t), where T c = diag(l, £,..., e r ~ 1 ). In these new 
variables, the first equation in (12) takes the form 

z-{t) = zf +1 (t), i = l,...,r- 1 

e~ r z°(t) = -kjx a (t) +v(t) 

ya(t) = z{{t). (17) 

Denoting rji(t ) = zf(t) — ai(t, x), i = 1, ... , r, where the stabilizing functions cti(t, x p ) are defined as 

= -Airyi(t) + di(£, x p ) - e T (t)Pb p g p (x p (t)) 

= -rji_ 2 (£) +d»_i(£, x p ), i = 3, . . . , r , 


OL2 (£, 3Cp) 
OLi (t , X p ^j 


( 18 ) 



and A,; > 0, * = 1, . . . , r — 1 are design parameters, the error system is written in the form 


eft) 


A m eft) + bp5p(£c p )77i(t) 


’ 0 ’ 


Vi(t) 

= 

— Ai? 7 i (t) +rj 2 ft) - e T (t)Pb p g p (x p ) 

+ 

0 

V rft) 

_ Tjr—l if) 


rjr—2(t) A r — i?/ r _i(t) 


_ 1 _ 



±(t) g 

£~ r f) r (t) = -kjx a (t) -£~ r a r (t,x p ) +v{t) (19) 

We define the control law as 

Vb(t) = kj X a (t) - X r T} r ft) - £~ r T]r-l(t) + £~ r Q. r (t, X p ) , (20) 

where A r > 0 is a design parameter. Then the closed-loop error system is reduced to 

x(t) = f{t,x(t))+gr] r (t) 

£~ r r)r(t) = -X r T] r (t) — £~ r iJ r -i(t) . (21) 

Lemma 1 The control signal (20) guarantees the boundedness of all closed-loop signals and the exponential 
convergence of the errors eft) and gift), i = 1, . . . , r to zero. 

Proof. Consider the following candidate Lyapunov function 

V 2 (e(t),r/(t)) =Vl ( e(t)) + r/?(t) + ■ ■ ■ + r/?(t) . (22) 

Its derivative can be readily computed to satisfy the inequality 

V 2 (t) = -e T (t)Qe(t) - Ai gift) X r -ig(._ 1 ft) - eX r g 2 r ft) < -gV 2 (t) (23) 

for all eft) G 1" and zft) G R, where g = A“* X (P) min [A“( n (Q), Ai, . . . , A r _i, -£ r X r ] with A max (P) denoting 
the maximum eigenvalue of the matrix P. Since V 2 ft) is radially unbounded, the inequality (23) implies 
the global exponential stability of the system (21). Since the reference input y c ft) is bounded, from the 
boundedness of eft) and x m ft) the boundedness of x(t) follows. Therefore, a(t,x) is bounded, implying the 
boundedness of z c ft). Therefore, the state x c a (t) is bounded. Then, the boundedness of internal state £ a (f) 
follows from the input to state stability of the actuator’s internal dynamics. □ 


Remark 1 We notice that the control signal Ubft) is comprised of two parts. The first part does not contain 
high gain and is designed from the perspective of retaining the slow mode of the actuator’s dynamics, since 
it cannot be altered without a high magnitude input signal. The second part is inversely proportional to the 
small parameter £ and is designed to cancel interconnecting terms. This part is a potential source of the high 
magnitude command that can exceed the actuator’s physical limits. 


Next, we apply singular perturbation method to the cascaded error system (19), and introduce a new 
(slow) time variable t = £ r t. Ignoring the actuator’s input to state stable internal dynamics and substituting 
r results in the following system 


®(t) 

dr 

dri r {T) 

dr 


f(r,x) +gvr(r) 


-kjxafr) +v(t) 


da r (r , x p ) 
dr 


(24) 


This is a singular perturbation problem. It is in the standard form if the algebraic equation 


f{T,x)+gi ]r (T) = 0 . 


(25) 


has isolated roots. Here we are concerned about the conditions of the existence of the isolated roots in 
general. For our purposes it is sufficient to note that the equation (25) has a root x* = [e = 0, r/i = 



0, ...,r ) r -2 = 0, r/ r -i = A r _iri r (t), which we denote by x* = h(t.,ri r ). Introducing a change of variables 
£(t) = x(t) — h(t, T] r ), the boundary layer system is obtained from the first equation in (24) as follows 

^ = (26) 
It is easy to see that the system (26) has an equilibrium at the origin. The reduced system takes the form 


%& = -*:x.(T)+»(r)-^2#^| (27) 

dr a aw W dr \e=h(T,rj r ) v > 

Unlike the previous case, here the control signal v(t) is designed to stabilize only the reduced system (27) 
and has the form 


U?(t) — k a X a (t . ) X r T] r (t) £ CX r [t., Xp) | e= h.(t,Ur) • (28) 

Applying the control signal v p (t) results in the following closed- loop reduced systems in the real time scale t 

Vr(t) = -e r Kr]r(t) , (29) 

which has a globally exponentially stable equilibrium at the origin (q r = 0). 

Lemma 2 Under the control action (28), the singular perturbation problem (24) with initial conditions 
eo, r/i (0) , . . . , Tj r (0) has a unique solution eft, e), ift(t, e), . . . , rj r (t, e), which satisfies the relationship 

rj r (t,e) - fj r (t) = O(e) (30) 

for all t > 0, where fj r (t) is the unique solution of the reduced system (27) with the initial condition i} r (o) = 
z(( 0) — a r (0, Xq). Further, there exists a time instance T* such that the relationship 

x(t, e) = 0(e ) . (31) 


holds for all t>T*. Moreover, all closed-loop signals are bounded. 

Proof. Consider the candidate Lyapunov function 

Vs(e(t),v(t)) = Vi(e(t))+rfi(t) + --- + rfi_ 1 (t). (32) 

Its derivative along the trajectories of boundary layer system (26) can be easily computed to satisfy the 
inequality 

V 3 (C{t),z{t)) < -pV 3 {C(t), z(t)) , (33) 

with /.t = A“) lx (P) min )A min (Q), Ai, . . . , A r _i] , which implies the global exponential stability of the boundary 
layer system. The smoothness and boundedness conditions of Theorem 1 follow from the assumptions 
imposed in this section. Therefore, from Theorem 1 it follows that (30) holds for all t > 0, and there exists 
a time instance T(e) such that the relationship 

x(t, e) — h(t, r] r (t)) = O(s) (34) 

holds for all t > T(e). From the exponential stability of (29) it follows that there exists a time instant T\ 
such that r] r (t) = 0(e) for t > T\. Then, from the definition of h(t,x(t)) it follows that the relationship 
h(t,x(t)) = 0(e) holds for t > Tf. Therefore (33) holds for t > T* , where T* = max[T(e), T{\. From the 
above relationship and the boundedness of the reference model (10) it follows that the plant’s state x p (t) is 
bounded, implying the boundedness of the stabilizing functions afit, x p ), i = 1, . . . , r. Since rj(t) is bounded, 
the boundedness of z c (t) follows, implying also the boundedness of the actuator’s state x c a (t). From the 
properties of the actuator’s internal dynamics it follows that the internal state £°(£) is bounded as well. □ 

Remark 2 Instead of the exponential tracking achieved by the backstepping control Vb(t) (20), the control 
law v p (t ) (28) achieves only £ -tracking for the original error dynamics, with a transient time greater than 
T*] . However, the latter control signal has a smaller magnitude than the former one. To see this, we first 
notice from the continuity of the function a(t , x) that a(t, x ) — a(t, *)| e =fc(t,z) = 0( £ )- Taking into account 
the relationship (80), we can write 

Vb(t ) - v p (t) = O(e) - £~ r rj r -i (t) (35) 

which implies that the two control signals differ by the high gain term £~ r ij r -i(t), which is inversely propor- 
tional to £ r and comes from Vb(t) (20). 



V. Control design for uncertain systems 


In this section, we solve the same control problem assuming uncertain plant’s and actuator’s dynamics 
that are also subject to external disturbances d p (t) and d a (t) respectively, which have bounded derivatives. 
To be precise, we assume that the parameter u>, matrices A and A a , vector b a , and continuously differentiable 
functions f p (x p ), g p (x p ), g a (y a ) are unknown. The control objective is to design a control law that uses the 
plant’s state x p (t) and the actuator’s output y a (t) for feedback to achieve tracking of the reference model 
(i°). 

As in conventional adaptive control, using the universal approximation theorem, we approximate the 
unknown functions by radial basis functions networks on compact sets fl x and Q ya 

f P (x p ) = wjif f (x p )+6f{x p ) 

g p {x p ) = w]ip g (x p )+6 g {xp) 

9a(y a ,Pa ) = wj<p a (y a ) +S a (y„) , (36) 

where Sf(x p ), S g (x p ) and S a (y a ) are the approximation errors, uniformly bounded by positive constants SJ, 
5* and 6* respectively. Using the relationships (36), the error dynamics (13) can be written as 

e(t) = A m e(t) + b p k T x p (t) + w]<p f (x p ) + wj ip g (x p )y a (t) + a p (t, x p , y a ) - vy c {t) , (37) 

where a p (t, x p , y a ) = S / (x p (t)) + S g (x p (t))y a (t) +d p (t) is the combined disturbance term, which is uniformly 
bounded by a positive constant a f as long as x p (t) £ and y a (t) £ 

In this section, we assume that the actuator’s dynamics has a relative degree m — r + 1 and is represented 
in the following observer canonical form 

s _1 a; Q (f) = A 0 x a (t) - qy a (t) +b a [u(t) + g a (y a (i),p a ) +d a (t)] 

Va(t) = cjx a (t), (38) 

where Aq and Cq are defined similar to definitions in the previous section with r replaced by ?n, q £ R m is an 
unknown constant vector, and b a = [0 ... b a ,m-r + 1 ■ • ■ ^om] T - Since the actuator’s dynamics is of relative 
degree m — r + 1, the b a ,m-r+i is not zero. Let it be positive. Taking into account the approximation in 
(36), the dynamics of x a (t) is written as 

e~ 1 x a (t) = A 0 x a (t) - qy a (t) + b a [u(t) +wjip a (y a ) +a a (t,y a )\ ■ (39) 

where cr a (t , y a ) = S a {y a ) + d a (t) is the combined disturbance term, which is uniformly bounded by a positive 
constant cr* as long as y a £ fl Va . 

To avoid generating high gain control signals we again apply the singular perturbation method. However, 
the method cannot be directly applied, since the exponential stability properties for the resulting boundary 
layer and reduced systems cannot be establish in the adaptive control framework without parameter conver- 
gence. The latter can be guaranteed only for sufficiently rich regressor or input signals, which in general is 
difficult to verify a priory. As an intermediate step, we follow the approach from Ref. 4 and design suitable 
prediction models for which the exponential stability properties can be guaranteed, and hence Tikhonov’s 
theorem can be applied. The control objective is met when the closeness of the prediction models to the 
corresponding dynamics is established. This will be done by the choice of adaptive laws independent of the 
control design. 

A. Prediction models 

We consider the following prediction model for the plant’s tracking error dynamics 

e(t) = A m e(t) + b k T (t)x(t) + wj(t)tp f (x) + wj (t)ip g (x)y a (t) + <r p (f)sign (i/j{t)) - vy c (t) , (40) 

where the initial conditions e(0) is chosen identical with e(0), = [e(t) — e(t)) T Pb p , and all variables 

with ”hat” are the estimates of the corresponding variables without ”hat”, and are available for the control 



design. Let the prediction errors be e(t) = e(t) — e(t). Its dynamics can be written in the form 


e(t) = A m e(t) + b p [k T + wj (t)ip f (x p ) +wj (t)ip g {x p )y a (t) 

+ (t)<P g (x P )y a (t) +v p {t,x v ,y a ) -<j p (£)sign(V>(£))] , (41) 

where we introduce parameter estimation errors k(t) = k — k(t ), Wf(t) = Wf — Wf(t), w g (t) = w g — w g (t), 

AW = A - AW- 

To be able to introduce a suitable prediction model for the actuator dynamics in (39), we re-parameterize 
them following conventional adaptive observer-estimator scheme (see for example 9,15 ) and introduce filters 

e _1 /3iW = (A 0 -b°kl)f3i{i) + b°y a {i), j 9 1 (0)=0 
£_1 /3 2 W = (A) - b°ko)(3 2 (t) + b°u(t), f3 2 { 0) = 0 
£ _1 /3 3 W = (A 0 -b 0 ko)/3 3 (t) + b 0 a a (t,y a ), /3 3 (0) = 0 

= (A-6°fc 0 T )S ( £)+6VIW, 3(0) = 0, (42) 

where b° = [0 ... 1 6 9 , t . . . 6 9 J T is chosen such that the pair (Aq, b°) is controllable and the polynomial 
s m-r + b° r+l s m ~^ + ■ ■ ■ + b° rn is Hurwitz, and the vector k 3 is chosen such that the matrix Aq — b°k 3 is 
Hurwitz. The actuator’s output y a (t) is represented as 

Va(t ) = #7^1 W + ^J^W + wju> 3 (t) + Of (t ) , (43) 

where £ K. m and = [b am . . . & a ,m-r+ i] T £ R r are unknown constant parameters to be estimated 
online, a>i(f) = /3i(£), w 2 (t) = [fai (t) ■■■ /3 2r (t)] T , ^3 (£) = cd"3(f) and Of{t) = cj (3 3 (t). Since the filters 
are stable from the boundedness of the inputs the boundedness of the outputs follow. That is Of(t) is 
bounded by some positive constant oj since o a (t,y a ) is bounded and the vector function u> 3 (t) is bounded 
since the RBFs are bounded. The prediction model for the actuator dynamics in (39) is introduced as 

AW = (t)uj 2 (t) + wJ(t)uj 3 (t) + o f (t)sign(ip(t )) , (44) 

where again the variables with ”hat” are the estimates of the corresponding variables without ” hat”), and 
AW = y a it) — AW is the prediction error, which is given by the equation 

AW = (t)vi(t) + W w 2 W + wj (t)u 3 (t) + o f (t) - Of (t)sign(^(f)) , (45) 

where A W = <&i — & 2 (t) = & 2 - AW and w a (t) = w a - w a (t). 

Since g(x p ) is assumed to be positive, the neural network approximation (36) can be chosen such that 
the unknown weights w g are positive. 11 The RBFs can be chosen to be positive as well. Then the product 
wj ( t)<p g (x p ) is always positive. Therefore, it is reasonable to design the adaptive law such that all compo- 
nents of w g (t) are positive. This will be done by means of the projection operator. 12 In the same fashion, 
since d 2 ,r > 0, the adaptive law for d 2 >r (f) is designed to keep it positive. Having this in mind, we design 
the adaptive laws for the estimates k(t),Wf(t),w g (t),o p (t) as follows 

AW = 7|V>WI 
k(t) = ^x p {t)ip{t) 

Wf(t) = 7 <p } (x)il>{t) 

Wg(t) = yH {w g (t), y a (t)ip g (x p )ip(t)} , (46) 

where 7 > 0 is the adaptation rate, n{-, •} denotes the projection operator, 12 introduced here to keep the 
estimates w g (t) bounded away from zero. 

The adaptive law for d 2 ,rW is defined by means of the projection operator. 

#2,rW = 7 n ($ 2 r(t), y a (t)u> 2 r(t) + ™ J (f )<p s (AV’ W w 2r (£)) . 


(47) 



The remaining adaptive laws are 


#1 (t) = 7 J/a(i)wi(t) +'ywj(t)<p g (x)ip{t)u 1 (t) 

$ 2 i(t) = 7 (y a (t)u 2 i(t) +wJ(t)ip g {x)ip(t)(j 2 i{t)j , * = 1, • • . ,r — 1 

^a(t) = 7y a (t)v3(t)+iw](t)<p g {x)i)(t)u 3 (t) 

Vf(t) = l\ya(t)\+^w](t)^p g {x)\%l}{t)\, (48) 

where the projection operator II (•, •) is introduced to keep u)(t) positive. 

The following lemma guarantees the closeness of the prediction models to the corresponding dynamics. 

Lemma 3 The adaptive laws (46), (4V and (48) guarantee boundedness of the estimates k(t),Wf(t),w g (t),a p (t) 
as well as w a (t), <r/(f) and the prediction errors e(t) and y a (t). Moreover, e(t) G £ 2 and 

Vait) G C 2 . 

Proof. Consider the following candidate Lyapunov function 

Vi{t) = e T (t)Pe(t) + y _1 [(< 7 * - a p (t)) 2 + fc T (t)k(t) + wj(t)w f (t) + wj(t)w g (t) 

+ 0(0!^) + 0(0 2 (t) +wj{t)w a {t) + (a* f -a fit)) 2 ] . (49) 

The derivative of V 4 (t) is computed along the trajectories of systems (41), and (46) and (48). 

V 4 (t) = e T (t){A(( % P + PA m )e(t) + 2 e 1 (t)Pb p [k T (t)x p (t) + w](t)<p f (x p ) + wj {t)ip g (x p )y a (t) 

+ (t)<fi g (x p )(dj {t)ui{t) + 0(t)u 2 (t) +wj(t)u: 3 (t) +<J f (t) - dy(i)sign(j/ 0 (t))] + a p (t,x p ,y a ) 

- o- p (t)sign(V>(i))] + 2y _1 [ - (a* - a p (t))a p (t) - 0 {t)k(t) - w J ( t)w f (t ) - wj (t)w g (t) 

- 0 ( 0 i (i) ~ ^ 2 ( 02 (t) - wj(t)w a (t) - (a} - a f (t))& f (t)] 

= -e T {t)Qe(t) + 2wJ (t)[i>(t)<p g (x p )y a (t) - 7 _1 &</(*)] + 2 wj (t)[i/>(t)tp f (x p ) - 7 _1 ™/ W] 

+ 2 0 (t)[ip(t)x p (t) — 7 _1 fe(t)] + 2[ip(t)a p (t, x p , y a ) - a p (t)\ip(t)\ - 7“Vp ~ a p (t))a p (t)] 

+ 20 (t)[if(t)w T g (t)ip g (x p )xi (t) ~ 7 _1 ^i (*)] + 2wJ (t)[i)(t)w](t)<p g (x p )u 3 (t) - r )~ 1 w a (t)] 

+ 2ti 2 (t)[ip(t)w] (t)(p g (x p )u> 2 (t) - 7 _1 ^ 2 (t)] + 2i/)(t)w](t)<p g (x p )a f {t) 

~ 2wJ (t)<p g (x p )a* f \if(t)\ + 2wJ {t)cp g (x p )[a* f - a f (t)\\ip(t)\ - 2y“ V/ - &f(t))a f (t )] . (50) 

Substituting the adaptive laws and taking into account the properties of the projection operator 12 

w](t) [y a (t)<p g (x p )il>(t) -n{w g (t), y a (t)ip g {x)ijj(t)}] < 0 (51) 

6 2 r (t) y a (t)u 2r (t) +w g (t)<p g (x)ip(t)cj 2r (t) - n|i? 2r (t), y a (t)u 2r (t) + w](p g (x)ip(t)u 2r (t)} <0, 

we obtain 

V 4 (t) < -e T (t)Qe(t) + 2ip(t)(Tp(t, x p , y a ) - 2a p (t)\'tp(t)\ - 2 y a (t) (i)uq(t) + (t)u 2 (t) + wj (t)u) 3 (t)] 

+ 2 )ip(t)w] (t)ip g (x p )a f (t) - 2wJ {t)<p g (x p )a}\ip(t)\ - ( a } - cr f (t))\y a (t)\ . (52) 

We notice that ip(t)a p (t,x p ,y a )-a* \ip(t)\ < 0, and iji(t)w0t)ip g (x p )af(t)-wJ(t)ip g (x p )a* f \'ip(t)\ < 0, since 
wj (t)tp g (xp) is positive by construction. Moreover, from the adaptive law in (48) it follows that a f(t) > 0 
as long as <j/(0) > 0. Therefore the following inequality holds 

(a) - d f (t)\y a (t)\ > a f (t)y a (t) - a f (t)sign(^(t))y a (t) . (53) 

Taking into account the ineciualities above and the error definition (45) we obtain 

V 4 (t) < -e T (t)Qe(t) -y 2 a (t), 


(54) 



which implies that the error signals e(f), k(t), ibf(t ), w g (t), a* — a p (t), (t), w a (t), 'd 2 (£) and a J — 

are bounded. Then, the estimates fe(£), Wf(t), w g (t), a p (t), $ 1 (f), w a (t), d 2 (t) and a f(t) are bounded as 
well. At this moment we cannot conclude anything about the boundedness of e(f), e(f), y a (t), y a {t ) or y a (t). 
However, we can integrate the inequality (54) and obtain 

[ [e T (r)Qe(r) + y 2 a (r)]dT < H 4 (0) - 14(f) . (55) 

Jo 

From Lemma 2 it follows that 14(f) is bounded, therefore e(f) £ £2 and y a (t) £ £ 2 - □ 


B. Adaptive control design 

We are interested in designing a control signal u(t) for the cascaded system 


e(t) = A m e{i) + b\k (t)x(t) + + Wg (t)<p g (x)y a (t) + & p (t)siga(ip(t)) - vy c {t) 

e "^(t) = (A 0 — b°kQ )f3 2 (t) + b°u(t ) , 


(56) 


where y a (t) = q T (t)x i(f) + b a (t)/3 2 (t) + tbj(£)% 0 (£) + <r/(f)sign(^(f)), such that e(f) — ■> 0 as f — > 00 . We 
notice that the second system in (56) is of relative degree r, stable and minimum phase. Therefore, it can 
be represented in the normal form 


e 1 dj 2 {t) = A 0 u> 2 (t) + b 0 [-kj (3 2 (t) + «(£)] 

e-'Cit) = A^(£) + b> 2 (£), (57) 


where the triplet (cj £ R r , Aq £ R rxr , 6 0 S R r ) is the normal form of representation of chain of r integrators, 
and the eigenvalues of AS £ (m-r) are ^g s t a ble ze ros of the minimum phase transfer function 

G(s) = Cq (ssl — A 0 + b°fe([)b 0 . The control design besides the tracking task must guarantee also the 
boundedness of u) 2 (t). Then the boundedness of £“(£) will follow from the input-to-state stability of the 
internal dynamics. 

To this end consider a change of variables z°(t) = T 0 w 2 (£), where T a = diag(l, £,..., £ r_1 ). In these new 
variables, the second equation in (56) takes the form 


m 


£ 


m 


z° +1 {t), i = 1, . . . ,r — 1 
~koT~ 1 z°{t) + u(t). 


(58) 


We notice that 


w] (t)<p (x)y a (t) = (x) q T + u(t)z°(t) + w a (f)x 3 (f) + d f (t)sign(tp{t)) 


where the inequality u>(t) > 0 is guaranteed by the adaptive law (48). Therefore, we design the stabilizing 
function z°(t) = di(f, x p ) for the tracking error prediction model as follows 

. u _ ^ _ -k T {t)x p (t)-w](t)<p f (x p )-a p (t)sign('ip(t)) + vy c (t) 

yt.X'p) — /s -r- 

d 2r (t)w] (t)ip g (x p ) 

q T (t)Xi(t) + wj{t)x 3 {t) + g/(£)sign(V>(t)) 

kr(t) 


If we introduce the variable r/i (f) as in the previous section, that is 771 (f) = z°(t) — ai(t,x p ), the error 
prediction model takes the convenient form 

e(£) = A m e(t) + bwj(t)(p g (xp)q 1 (t) . (60) 

However, the dynamics of 771 (f) involves the derivative of &i(t, x p ), which is not available in this case because 
of unknown terms in the plant’s dynamics. One remedy for the situation is the command filtered approach 



presented in. 3 Following Ref., 3 we first filter the stabilizing function x p ) through a second order stable 
filter 


zn{t) = u> a z 12 {t) 

zi 2 (t) = - 2 C, a L 0 a zi 2 {t) - u a [zu(t) - ai(t, x p )\ , (61) 

then use the error variable = z°(t) — z\\(t) instead of iji (t). The tracking error prediction model now 
can be written as 


e(t) = A m e(t ) + b p wj (t)<p g (x p )f)i(t) + b p wj (t)ip g (x p )[z n (t) - an(t, x p )] , (62) 

The compensated error e(t) = e(t) — £(t), where £(t) is defined according to equation 

€(t) = A m £(t) + b p wj (t)<p g (xp)[z n (t) - an (t, x p ) + £ a i(t)] , (63) 

satisfies the equation 

Ht) = A m e(t) + bpwj (t)<p g (x p )fji (t) , (64) 

where fji{t) = fji(t) — tai(t) is to be defined in the next step. 

Following the block-backstepping' and command filtering 3 procedures, the stabilizing functions d, ,* = 
2, . . . , r are defined as follows. 

a 2 (t,x p ) = — Ait)i ( t) + u a z 12 (t) - (e(f) - £{t,)) T Pb p wJ (t)ip g {x p ) 

Ui{t,x p ) = - fji~ 2 {t) + u> a Zi-i, 2 (t ) , « = 3, . . . , r , (65) 

where ?),:(t) = z°(t) — zu(t), za(t) is generated through the stable command filter 

Zn(t) = L0 a Z i2 (t) 

= -2CaU a z i2 (t) - u> a [zn(t) - ati(t,Xp )\ . (66) 


The prediction models (56) in terms of error signals , fj T (t) can be written in the following compact 

form 


1 

1 


A m e(t) + bpibj (t)<p g (x p )fji (t) + b p w] (t)<p g (x p )[zu(t) -&i(t,x p )] 


’ 0 " 



= 

-Xifji - (e(t) - £(t)) T Pb p wJ {t)<p g (x p ) + Z 2 i(t) - a 2 (t,x p ) 

+ 

0 

Vr(t) 

_ rjr-l(t) 


-T) r - 2 (t) - X r -ifj r -l(t) + z r i(t) - a r (t,x p ) 


_ 1 _ 



x(t) 9 

£~ r fjr{t) = -s~ r z rl (t) +u(t) - ko T~ 1 z°(t) (67) 


The compensated error is defined as r/i(t) = r)i(t) — £ a i{t), where £ a *(f) is generated by the dynamics 


^az(^) — ^iCai(^) T ^i+l,l(t) Qfj(t, 3?p) T ^a,i+l (^)j ^ — 1 5 ^ 1 

L,r = 0. (68) 

From the above constructions it follows that fj r (t) = rj r (t). The dynamics of the compensated error can be 
written as 


= — A1T71 (t) + 772(f) - e T (t)Pb p w]{t)ip g {x p ) 

= - Xifji(t) +fj i+ i(t), i = 2, . . . ,r — 1 

= —£~ r z r \{t) + u(t) — k$ T~ 1 z°(t) . 


m (t) 

£~ r fjr(t) 


(69) 



We apply the singular perturbation method to the cascaded system 


e(t) 


A m e(t) + b p wj (t)<p g (x p )r} 1 (t) 


’ 0 ’ 



= 

-Aim (t) + m{t) - e T {t.)Pb p w] {t,)ip g {x p ) 

+ 

0 

Vr(t) 



-fj r - 2 (t) - Xr-lfjr-l{t) 


_ 1 _ 



s(t) /(t,x) a 

£~ r fjr(t) = —E~ r Z r i(t) + u(t) — fej" T~ 1 Z°(t) (70) 

by introducing slow time variable r = e r t. In the slow time the systems are written as 

= +gf] r (t) 

= —+u(t)-koT~ 1 z 0 (t) (71) 

This system is in the form of the standard singular perturbation model, if the algebraic equation 

f(t,x)+gf} r (t) = 0 (72) 

has isolated roots. We are not concerned about finding all the roots of the equation (72), but we are interested 
in the root x* = [e = 0 fji = 0 ... fj r -2 = 0 fj r - i = A“_ 1 r) r (t)] T , which we denote by x* = h{t,fj r ). 

Introducing a change of variables C( T ) = *(t) — h(t,fj r ), the boundary layer system is obtained from the 
first equation in (71) as follows 

= f(t,C + h(t,fj r ) + gfj r (t) = f{t,C) ■ (73) 

clt 

The system (73) has its equilibrium at the origin. The following lemma shows that it is globally exponentially 
stable. 

Lemma 4 The boundary layer system (73) is globally exponentially stable. 

Proof. Consider a candidate Lyapunov function 

V 5 (t) = e T (t)Pe(t) +fjl(t) + -.-+ f r _, (f) . (74) 

It is straightforward to compute its derivative along the solutions of system (73). After some algebra we 
obtain 


r dx(r) 
dr 

dfj r (t) 

dr 


V 5 (t) = -e T {t)Qe{t) - Au 7 i(i) K-iVr-i(t) ■ (75) 

Denoting fi = A“) lx (P) min [A min (Q), Ai, . . . , A r _i] the following inequality can be readily derived 


V 5 (t) < -fj,v 5 (t ) , (76) 

which along with the radial unboundedness of V$(t) implies the global exponential stability of (73). □ 

The reduced model is derived from the second equation in (71) by substituting x (r) = h(r, fj r ) and s = 0. 
Since the only term fej T“ 1 z°(t) does not explicitly depend on x , we can write the reduced system in the 
form 


dijrjr) 

dr 


dz r i ( t ) 
dr 


+ v(t)-k( j T 0 1 z°(t). 


(77) 


The control signal v(t) is designed in the actual time scale as 

u(t.) = -A r fjr{t) + k^T~ l Z 0 {t) + £~ r U> a Z r2 {t) , 


(78) 


which in the actual time scale results in the following closed-loop reduced system 

Vr(t) = -£ r Kv r(r). 


(79) 



Obviously the reduced system has a globally exponentially stable equilibrium at the origin. The rest of the 
conditions of Theorem 1 follow from the smoothness and boundedness assumptions imposed on the reference 
input, as well as on the plant’s and actuators’s dynamics. Therefore, for allt > 0 the singular perturbation 
problem (71) has a unique solution x(t,e), fj r (t,e) with initial conditions ai(0) and fj r ( 0)) respectively, and 
the relationship 


T) r {t,e) - rj r (t) = 0(e) (80) 

holds, where is the unique solution of the reduced system with initial condition fj r ( 0) = z°( 0) — z r i(0). 
Moreover, there exists a time instance T(e) such that the relationship 

x(t, e) — h(t, fj r {t)) = O(e) (81) 

holds for all t > T(e). From the exponential stability of (79) it follows that there exists time instant Tf such 
that fj r (t) = 0(e) for t > T\. Then, from the definition of h(t,fj r (t)) it can be seen that the relationship 
h(t,fj r (t)) = 0(e) holds for t>T\. Therefore, 

x(t,e) = 0(e), (82) 

holds for t > max[T(e), Tj]. Next, we apply the result of Ref. 3 to the error systems (19) and (67). Since 
x(t,e) = 0(e) and fj r (t) = O(e) for all t, > max[T(e), Ti], it can be shown that 

x(t,s) = O(e) + Oiu- 1 ) 

Vr{t, e) = 0(e) + 0(ix~ 1 ) (83) 

for all t > max[T(e), Ti], where [x(t,e) fj r (t,e)] is the solution of the error prediction system (67). 

From (83) it follows that e(t) and f]{t) are bounded. Then, from Lemma 2 it follows that e(t) is bounded. 
Since the reference command y c (t ) is bounded, x m (t) is bounded, implying the boundedness of x p (t). Since 
Lemma 2 guarantees the boundedness of the parameter estimates, it follows that the stabilizing functions 
cq(t, x p ), i = 1, . . . , r are bounded, implying also boundedness of (zn(t), Zait)), i = 1, . . . ,r. Then from the 
definition of error signal fj{t) it follows that z°(t ) is bounded, implying the boundedness of y a (t ) and /3 2 (0- 
Then the control signal u(t) is bounded as well. From the the definition of y a (t) (44) it follows that (t) is 
bounded, hence (3 1 (t) is bounded. Then from the filter equation (42) it follows that y a {t ) is bounded. That 
is, all closed loop signals are bounded. Therefore, e(t) and y a {t ) are bounded. From the Barabalat’s 14 it 
follows that e(t) — > 0 and y a (t ) — > 0 as t — > oo. 

The obtained results are formulated in the form of the following theorem. 

Theorem 2 Consider the uncertain plant (8) and slow actuator (9). The adaptive controller given by (28), 
along with the stabilizing functions (If) and (18), auxiliary filters (61) and (66), the prediction models (40) 
and (44), and the adaptive laws (46), (41) and (48), guarantees that the plant’s output tracks the output of 
the given reference model (10) with the error directly proportional to the small parameter in the actuator’s 
dynamics and inversely proportional to the frequency of the filters (61). 

VI. Simulation Example 

We demonstrate the benefits of the approach on a first order unstable system 

x p {t) = Xp(t) +x a (t), (84) 

controlled by the actuator 

x a (t) = -0.1 (x a (t) - u(t)) . (85) 

For this simulation, the control objective is to track a reference command y c (t) = 1. That is e(t) = x p (t) — 1. 
The control signal Ub(t) takes the form 

= x a {t) + a m e(t) + x 2 p {t) 

= ~ kaV(t ) + x a (t) - e _1 e(t) - e _1 { a^e(f) - a m r](t ) - 2 x p (t) [x 2 p {t) + x a {t)] } , 


v(t) 

Ub(t) 


( 86 ) 



Plant output 



Figure 1. Output tracking performance for different control strategies. 


whereas the signal u p (t) has form 

Up{t) = -k a ii(t) + x a (t ) + £ -1 <2 (— r/(t) + l) + X a (t) ( ——r](t) + 1 

Figure 1 displays the tracking performance for three controllers: feedback linearization, backstepping 
design and proposed design with time scale reparation. For the backstepping design and proposed design 
the control gains are the same: a m = 3, k a = 1. For the feedback linearization the control gains are chosen 
such that the eigenvalues of the resulting error systems are the same. As it can be seen from the plots the 
tracking is almost the same. The actuator’s output performance is displayed in Fig. 2 and the control signals 



Actuator output 



Figure 2. Actuator’s output performance for different control strategies. 

are displayed in Fig. 3. As it can be seen the magnitude of the control signal generated by the singular 
perturbations method based approach is substantially smaller the that of generated by feedback linearization 
and backstepping design. 

We notice that with the control signal w p (f), the error system that is comprised of the boundary layer 
system and reduced system has eigenvalues —3 and —0.1. Next we run the simulation for the feedback 
linearization and backstepping approach with the gains chosen such that the error system has poles at the 


Control signal 



Figure 3. Control signals for different control strategies. 


same location. The tracking performance and the corresponding control signals are displayed in Fig. and 
Fig. respectively. As it can be seen from the figures the tracking is sluggish for both designs, as it is expected, 
though the control signals have smaller magnitude. 
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Figure 4. Output tracking performance for backstepping and feedback linearization with poles at —3 and —0.1. 


VII. Conclusion 

We considered a problem of controlling a plant with a slow actuator. It has been shown that the 
singular perturbations method can be applied for this problem. The resulting control signal has substantially 
smaller magnitude than one generated by conventional design technique, although the tracking error can be 
guaranteed to be only bounded with the bound proportional to the small parameter in actuator’s dynamics. 


VIII. Acknowledgment 

The authors wish to thank Dr. Eugene Lavretsky of The Boeing Company for fruitful discussions. 


3 


Control signal 


2.5 - 

2 f 

J 

1.5 \ 

1 V 

\ 

0.5 - ' 

\ 

\ 

0 - V 

Y 

- 0.5 - 


-1 L 
0 


- - - Feedback linearization 
Backstepping 


5 


10 


15 


Figure 5. Control signals for backstepping and feedback linearization with poles at —3 and —0.1. 


References 

1 C.-C. Chen. Global Exponential Stabilization for Nonlinear Singularly Perturbed Systems. IEE Ptoc. - Control Theory 
Appl. , 45(4):377-382, July 1998. 

2 P. D. Christofides and A. R. Teel. Singular Perturbations and Input-to-State Stability. IEEE Trans. Autom. Contr ., 
41 (11): 1645-1650, Nov 1996. 

3 J. A. Farrell, M. Polycarpou, M. Sharma, and W. Dong. Command Filtered Backstepping. IEEE Trans. Autom. Contr ., 
54(6):1391-1395, June 2009. 

4 N. Hovakimyan, E. Lavretsky, and C. Cao. Adaptive Dynamic Inversion via Time-Scale Seperation. IEEE Trans. Neural 
Networks , 19(10):1702-1711, 2008. 

5 H.K. Khalil. Nonlinear Systems, Third Edition. Prentice Hall, New Jersey, 2002. 

6 P. V. Kokotovic, H. K. Khalil, and L. O’Reilly. Singular Perturbations Methods in Control: Analysis and Design. 
Academic Press, New York, 1986. 

7 M. Krstic, I. Kanellakopoulos, and P. Kokotovic. Nonlinear and Adaptive Control Design. John Wiley & Sons, New 
York, 1995. 

8 C. L. Lin and B. S. Chen. On the Design of Stabilizing controllers for Singularly Perturbed Systems. IEEE Trans. 
Autom. Contr., 37(11):1828-1834, 1992. 

9 K.S. Narendra and A.M. Annaswamy. Stable Adaptive Control. Prentice Hall, 1989. 

10 p. Sannuti and A. Saberi. Special Coordinate Basis for Multivariable Linear Systems - Finite and Infinite Zero Structure, 
Squaring Down and Decoupling. The International Journal of Control, 45:1655-1704, 1987. 

11 J. Park and I. Sandberg. Universal Approximation Using Radial Basis Function Networks. Neural Computation, 3:246- 
257, 1991. 

12 J. B. Pomet and L. Praly. Adaptive Nonlinear Regulation: Estimation from the Lyapunov Equation. IEEE Trans. 
Autom. Contr., 37(6): 729-740, 1992. 

13 A. Saberi and H. Khalil. Qudratic-type Lyapunov Functions for Singularly Perturbed Systems. IEEE Trans. Autom. 
Contr., AC-29:542-550, 1984. 

14 S. S. Sastry and M. Bodson. Adaptive Control: Stability, Convergence and Robustness. Prentice Hall, 1989. 

15 J.J. Slotine and W. Li. Applied Nonlinear Control. Prentice Hall, New Jersey, 1991. 

16 A. N. Tikhonov. On the Dependence of the Solutions of Differential Equations on a Small Parameter. Mat. Sborni. N. 
S. 22(64) (i n Russian), pages 193-204, 1948. 


